Spin fractionalization of an even number of electrons in a quantum dot 
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Kondo resonant conductance can occur in a quantum dot with an even number of electrons A'' in 
the Coulomb blockade regime, with bottom and top contacts attached to it in a pillar configuration, 
in a magnetic field B_i_ = B* along the axis. The field B* is tuned at the crossing of levels with 
different total spins and the spin degeneracy of the dot energies should be lifted. The case of the 
degeneracy of the singlet-triplet states in a dot with A'^ = 2 is analyzed in detail. Coupling to the 
contacts is antiferromagnetic due to a spin selection rule, so that total spin on the dot is 5 = 1/2 
in the Kondo state. 
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Quantum dots {QD) are remarkable because of 
Coulomb oscillations R]. At very low temperature, lin- 
ear conductance (at vanishing bias Vsd) is zero, except 
for peaks at discrete values of the gate voltage Vg, when 
it is energetically favourable to add one extra electron to 
the dot. Therefore Vg controls the number of particles 
on the dot in the Coulomb Blockade (CB) regime: its 
charge is quantized and so is the spin. Dots with an even 
(odd) number of electrons have integer ( half integer ) 
total spin. 

However, by lowering the temperature further below 
some value Tr-, finite conductance can be found exper- 
imentally at zero voltage [^-Q, relatively insensitive to 
small changes of Vg within the CB regime, as a meso- 
scopic realization of the physics of the Kondo effect [^Q . 
The zero voltage anomaly in the conductance is only seen 
when population is odd [gj. While the charge remains 
quantized, the spin quantization of the correlated state 
between the dot and its contacts is changed, due to the 
formation of a Kondo singlet involving the odd electron. 

A weak magnetic field parallel to the dot plane, i3||, 
splits the peaks in the differential conductance into two, 
located at Vsd = ±5*Mb-S||i (where g* is the effective 
giromagnetic factor for electrons in this geometry and 
fiB = eh/2mc is the Bohr magneton), until eventually 
they disappear ^^■ 

On the contrary, a magnetic field B±, orthogonal to 
the dot area, is believed to be much more destructive, 
due to orbital changes of the electron states |§] , in addi- 
tion to the lifting of spin degeneracy , that prevents the 
formation of the Kondo correlated state. 

We suggest here the possibility that, due to accidental 
crossing of the lowest non degenerate many body energy 
levels of the QD at special values of B± = B* , orthogonal 
to the dot connected to the leads in a pillar structure, the 
Kondo resonance can still be present, irrespective of the 
parity of the electron number N p,pi3|. In particular 
we discuss here the case of A^ = 2 and show that under 
suitable conditions a correlated state sets in, which is 



characterized by fractionalization of the total spin. 

Indeed a QD at high _Bj_ (in the z — direction) displays 
a rich level structure that can be monitored attaching 
bottom and top contacts to it in a pillar structure via 
linear and non linear transport |ll|-[l4|]. Assuming the 
QD to be effectively two-dimensional and confined by a 
parabolic potential, its states in isolation can be labelled 
by the total spin S and the z— component of the total 
angular momentum M and spin Sz- In increasing Bj_, 
many level crossings can occur for an N particle QD 
between states of increasing Af , as well as states with 
higher total spin S. 

It is well known that the ground state (GS) of an iso- 
lated dot with N = 2 electrons changes from singlet (^S' : 
5 = 0, M = 0) to triplet ^Ts^ : S* = 1, M = 1) at a value 
B±^ = B* , due to crossing of levels ||l^ . Correspondingly, 
exact diagonalization shows that the A^ = 3 electron dot 
is in a doublet state {^Ds^ : S = 1/2, M = 1 or 2 de- 
pending on the field value) (see fig.0) Q. As shown 
quantitatively in the following, the dot parameters can 
be chosen in such a way that B± is high enough to lift 
spin degeneracy. We take Vg tuned in the CB region be- 
tween N = 2 and A^ = 3 electrons on the dot, and limit 
our discussion to four states : 1) dot with A^ = 1 in ^D, 
2-3) dot with A^ = 2 in one of the two degenerate states 
2 5*0, and ^Ti, 4) dot with A^ = 3 in ^I^i/s. Taking only 
these states allows the mapping of the problem onto the 
non degenerate Anderson model with the following levels: 
1) empty level of zero energy, 2-3) singly occupied level of 
energy Cd (doubly degenerate), 4) doubly occupied level 
of energy 2ed + U. Here U is very large (see below) and, if 
levels 2-3 and 4 are symmetric with respect to the chem- 
ical potential of the contacts fi, what can be controlled 
by Vg, average occupancy of levels 2-3 is n^ = 1. This 
corresponds to A^ = 2 electrons on the dot. 

The key observation is the following spin selection rule 
for the isolated dot [Q : because the A^ = 3 GS for 
the dot only has Sz = 1/2, hybridization of the dot with 
the contacts can only occur by virtual occupancy of the 



dot with an extra spin up electron when the dot is in 
the N = 2 singlet state ^S" or with an extra spin down 
electron when the dot is in the N — 2 triplet state with 
Sz = I {^Ti state). Although spin degeneracy of the dot 
is assumed to be lifted, the physics is the same as that 
of the single channel Kondo model with effective spin 
S'^-^^ = 1/2, in the absence of magnetic field. Here is 
S^-^^ — Sz — 1/2 and coupling with the delocalized elec- 
trons is antiferromagnetic due to the spin selection rule 
mentioned above. Provided hybridization to the contacts 
A is large enough and temperature is below the Kondo 
temperature Tk, the system flows to the strongly corre- 
lated fixed point. Formation of the Kondo singlet corre- 
sponds in this case to an average spin on the dot frozen 
at 5*2 = 5, that is to fractionalization of the total spin of 
the dot, together with a strong zero bias anomaly in the 
conductance. 

In this letter we discuss the properties of the GS, bor- 
rowing from the known results of the non degenerate An- 
derson model [^. In particular, we discuss the spin frac- 
tionalization, the practical realizability of the system and 
the robustness of the GS versus small changes of the mag- 
netic field B± with respect to the value B* that makes 
the states IS") and \T) degenerate in energy. 



tacts within the CB valley of the conductance at A^ = 2 
(/iAT < /I < iiN+i = ^'^'^E —^Ei). Electrons in the con- 
tacts arc non interacting and have an energy dispersion 
Cka involving just one single Landau level and linearized 
around /i with Fermi velocity vp„. The wavevector k is 
in the z— direction and the two contacts are assumed 
to be equal (Hl/a = Y.ka «^fe,<^^L/fi kcr^L/R ka) and three 
dimensional enough, in order to have both spin polariza- 
tions at /i. 

The dot Hamiltonian Hd, restricted to the Fock space 
with A^ = 1,2,3 electrons can be represented in terms of 
the occupation numbers rij for the degenerate levels ^-B^: 
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FIG. 1. Calculated energies of a circular quantum dot vs 
magnetic field Bx (expressed in cyclotron energies fiojc) for 
N = 2 (a),(c) and N = 3 (b),(d). Here the level spacing of 
the confining potential is hujd = 4meV and the strength of the 
Coulomb interaction is Ec = SmeV (notations can be found 
in the text). (c),(d) magnify the marked area of (a),(b), 
respectively 



We denote the two dot degenerate energies by '^Ei 
{i — I : triplet state, i = 2 : singlet state) and take 
the chemical potential /i of the left and right bulk con- 



Hd = M2 ^ ni + Unin2 



(1) 



where U — 113 — 2/i2 is the charging energy for adding the 
third electron to the dot. Here is n^ = clci, the vacuum 
of the c operators is the dot state with one electron and 
the c^ create the two degenerate states by acting on it. 

The model Hamiltonian is H = Hl + Hr -\- Hd -\- Ht, 
where Ht is the hybridization term acting only at z — 0. 
According to the spin selection rule stated above, it takes 
the form: 
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The tunneling amplitudes F^/flj-cri, which depend on 
the actual features of the device, have been taken all 
equal for simplicity. 

We start integrating out the contacts fields (fi, = 1 in 
the following): 



Z(m)« fu,(^Dc,Dcl^ 



-Aa „-l3\r\^-£^^,,clK{iuj^)ci 



(2) 



with 



Ad 



^He 



f (J , \ t 



\...j. 



Uc[ciC2C2 



and K{ujm) = X]fc(*'^'" + ^^) ^ ^^ taken independent 
of a for simplicity. In the limit of large U we have: 



exp[/ dT[Unin2 ~ fj.{ni + n2)]] 
Jo 



= e^'^' S{ni + 712 - 2fi/U) ■ e"* ^0 drim-n^f ^ 

where the delta function implements the constraint of 
single site occupancy in the symmetric case, fl — ji — ^12 = 
U/2. 

We introduce X{t) as an auxiliary boson field coupled 
to the differential occupancy of the two degenerate states 
by means of an Hubbard- Stratonovitch decoupling: 



Then, the partition function of eq.(S)reads: 
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where is G,l{iujm) = itj„i + |rp_?ir(icj„j). We will 
take the saddle point solution of eqj^. At saddle point 
(see below) the constraint is satisfied in the average over 
the time (3 because {n, - \) ~ ~{-iyXs.p.{T)/2U = 
and the integrals for i — 1,2 decouple, giving a contri- 
bution to the effective action for the field X that, for 
1, is given by J^ dg J^ dT[G(gX),i{T,T+) - ^]X{t) 
G^o5 - gX i X ^ -X ior 1 ^ 2) 
|_8[. This representation allows to recognize the dou- 
bly degenerate GS of an effective spin S'^^f = 1/2 with 
Sff = {S^ - 1/2), driven by the field X{t) and pro- 
duced by virtual tunneling of electrons on and off the 
dot at energy /i. 

We now follow ref. [Q quite closely, to map eq.|| on 
a 1-dimensional Coulomb gas (CG). We define ^(r) = 
X{t)/A, where A = 27rA^(0)r2 and N{0) = L/ttvf {L is 
the size of each single contact) is the density of states at 
the Fermi energy. At low temperatures, saddle point so- 
lutions of the resulting single particle effective action are 
sequences of instantons ^±{t,1) = ±^otanh((r — ti)/To) 
(ti are the centers, / = 1, 2, . . .) corresponding to jumps 
between the two minima of the effective potential, V[£] — 
A^/U^^ - 2A/7r[^tg-i(0 - 5ln(l + ?^) |l|, located at 
^0 = ±t//2A, that interact via a logarithmic potential 
a^ln\{U~t,)/To\. 

The bare strength of the logarithmic interaction is 
= 2(1 - AA/Utt^). The bare fugacity of the CG is 

9[ lb = To cxp(— y^) {A ^ tqU is the action of one 
single blip). The scaling of the fugacity and the renor- 
malization of the coupling constant induced by processes 
of fusion of charges lead to the Renormalization Group 
(RG) equations |2(][| : 
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The flow is towards F — > cx) and a^ 
the scaling invariant energy is Tk 
{UA)2e~^ c//(4A)^ Condensation of instantons in the 
doubly degenerate GS leads to the Kondo singlet < 
5"=// >= or < 5 >= i on the dot. 



An heuristic argument for the fractionalization of the 
total spin S of the dot reads as follows. 

The dynamics of the field ^(t) of action A can be mim- 
icked by a two level system hamiltonian H21 with hop- 



ping energy A 
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The role of the interaction is to project out higher en- 
ergy components from the dot states. Denoting by |± > 
the two eigenstates of H21, the dynamics of the dot is 
between states U9\s > /7V|(.9) and U3\T > /J\f^{g), 
where \S > and \T > are the states localized in the 
two potential wells and J^s/xig) ^-re their normalizations, 
n^ = 1 — gll^ is the Gutzwiller operator partly eliminat- 
ing the component on the high energy eigenstate |— > 

(.9e(0,i))@. 

Now we state the connection with the CG picture given 
above. The frequency A/27r is obviously related to the 
average number of flips: /3A/27r =< N >= Ftanhy 
and probability of having the system in states |± > at 
temperature l3~^ are 
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(5) 



Equating these probabilities to the ones given by P± — 
|< ± |n9| S'/T >|^/7V5/t(5) yields the correspondence 

1 _ ^ ^ -2-!tY tanhY 

The z— component of the total spin of the projected 
singlet and triplet state are given by: 
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Because Y scales to inflnity, 5 — > 1: the higher energy 
state completely decouples and < Sz >s/T,g^ \- 

We now give some estimates for a possible practical 
realization. 

If the level separation of the parabolic conflning po- 
tential of the dot is Twod = AmeV , the effective fre- 
quency in presence of Bj, is ujq = \/oj^ + ^c/^ (^c — 
eB±/m*c is the cyclotron frequency with m* = 0.067?7ie 
for GaAs). The Coulomb interaction between electrons 
in the plane is Ec x lo/\r -- r'l ( where Ec — 3meV and 
Iq — h/m*ujo). This choice of the parameters gives a 
singlet-triplet (S-T) transition for the dot with N = 2 
at B± = _B* sa 4.5 Tesla, close to one of the Delft 
experiments |ll|. Fig.^ shows the total energy of the 
dot vs magnetic field for A^ = 2 and A^ = 3. Zeeman 
spin splitting of the energies is not included in the pic- 
ture for clarity, although it is sizeable at these fields: 
g*liBB = O.llrney.Here we assume the bulk giromag- 
netic factor for GaAs g* = —0.44, but it could be lower 
P2| , P3| . In any case, at temperatures below 300m AT the 
spin degeneracy is lifted. According to fig.^ the energy 
difference between the GS's when three and two elec- 
trons are on the dot is U ='^ E —^ E — 12.5meV at the 
S-T crossing, which is about 1/3 higher than the value 



Ec + hujo — i.meV + S.SSmeV predicted by a constant 
interaction model. The TV = 3 particle state is a doublet 
{S = 1/2 ) ^D with M = 1. Higher magnetic fields pro- 
duce its crossing with the doublet with Af = 2, while, in- 
creasing B_\_ further, the quadruplet state {Q : S = 3/2) 
becomes lower in energy. 

The order of magnitude of the hybridization of the 
dot with the contacts A is fixed by the requirement 
that neighbouring CB peaks in the conductance do not 
overlap 0. Their separation is ~ hajQ. For values of 
A ^ (0.1 — 0.2)huJo, Tk (given below) is of tenths of 
mK . We have neglected the excited states of the two 
particle dot in the description. This is allowed because 
Tk is much less than the spin splitting of energy levels. 
The next available state in the spectrum of the isolated 
dot with iV = 2 is the triplet state with S'z = ( ^Tq). 
Because it is non degenerate, no Kondo-like coupling is 
possible. This state can only hybridize with conduction 
electrons via cotunncling processes which, in the elastic 
case, give rise to currents that are linear in the biasing 
voltage Vsd and of the order of ^ (A^/[/) Vsd- 

We conclude that the problem is equivalent to a non 
degenerate Anderson model in absence of magnetic field. 
At B± = B* , the N ^ 2 particle dot can be in the "^Sq 
or ^Ti state, but addition of one more electron requires 
an extra energy U. Furthermore is U >> A. The results 
are expected to be quite robust with respect to small 
changes of the magnetic field away from B* , especially 
for B± < B* because the reference field i^o — (C^A) 2 //i^ 
is of the order of Tesla. For Bj^ > B*, instead, the state 
^Ti will be eventually favoured in energy. 

A small a.c. magnetic field |i3||| << \5B\ ^ Tk/jJ-b 
( with 5B = B± - B* < 0) could induce transi- 
tions between the two degenerate Kondo GS 's, with 
Rabi frequency hujR = Md\SB\. The quantity Md/2 « 
29* f^B [1 + SB/{ttTk)] is the expected magnetic moment 
of the Dot with A^ = 2, in the Kondo state. The measure- 
ment, if feasible, should be performed just off B± — B* 
for SB < 0, not to match the energy separation of the 
two levels ^Ti and ^Tq. 

In conclusion, we have discussed the possibility that 
the total spin S of a dot in CB with N=2 is strongly cou- 
pled to the spins of the electrons in the leads. This kind 
of mesoscopic Kondo effect would give rise to fractional- 
ization of S in the dot, i.e. to a sort of 'spinon box'. 
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